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Abstract 

We review the recent literature on lattice simulations for few- and many-body systems. We 
focus on methods that combine the framework of effective field theory with computational lattice 
methods. Lattice effective field theory is discussed for cold atoms as well as low-energy nucleons 
with and without pions. A number of different lattice formulations and computational algorithms 
are considered, and an effort is made to show common themes in studies of cold atoms and low- 
energy nuclear physics as well as common themes in work by different collaborations. 
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I. INTRODUCTION 

In this article we review the literature on lattice simulations for few- and many-body 
systems. We discuss methods which combine effective field theory with lattice methods and 
which can be applied to both cold atomic systems and low-energy nuclear physics. Several 
recent reviews have already been written describing quantum Monte Carlo methods for a 
range of topics. These include Monte Carlo calculations in continuous space for electronic 
orbitals in chemistry [l|, solid state materials superfluid helium and few-nucleon 

n 

systems |4J . There are also reviews of Monte Carlo lattice methods for strongly-correlated 
lattice models 5], lattice quantum chromodynamics at nonzero density [6[], and a general 
introduction to lattice quantum chromodynamics [TJ. 

Lattice simulations of quantum chromodynamics (QCD) are now able to accurately de- 
scribe the properties of many isolated hadrons. In addition to isolated hadrons, it is 



3 



also possible to calcu 



ing 



ate low-energy hadronic interactions such as meson-meson scatter- 



131 ] . Other interactions such as baryon-baryon scattering are com- 



putationally more difficult, but there has been promising work in this direction as well 



Q 



15 



13, 



18, 3. 



201 ] . A recent review of hadronic interaction results computed from 



lattice QCD can be found in Ref. 211 ]. 

However for few- and many-body systems beyond two nucleons, lattice QCD simulations 
are presently out of reach. Such simulations require pion masses at or near the physical 
mass and lattices several times longer in each dimension than used in current simulations. 
Another significant computational challenge is to overcome the exponentially small signal- 
to-noise ratio for simulations at large quark number. For few- and many-body systems in 
low-energy nuclear physics one can make further progress by working directly with hadronic 
degrees of freedom. 

There are several choices one can make for the nuclear forces and the calculational method 
used to describe interacting low-energy protons and neutrons. For systems with four or fewer 
nucleons, a semi-analytic approach is provided by the Faddeev-Yakubovsky integral equa- 
tions. Using this method, one study [22| looked at three- and four-nucleon systems using 



the Nijmegen potentials 



23 



with the Tucson-Melbourne 



, CD-Bonn potential 



2J], and AV18 potential [25(, together 



261 ] and Urbana-IX 27j three-nucleon forces. A different inves- 



tigation considered the^same observables using a two-nucleon potential derived from chiral 

Another recent study 29[ considered the low-momentum interac- 



effective field theory 
tion potential Vj ow k 



J29- 

2J31 



3l| . This method used the renormalization group to derive effective 
interactions equivalent to potential models but at low cutoff momentum. 

For systems with more nucleons approaches such as Monte Carlo simulations or basis- 
truncated eigenvector methods are needed. There is considerable literature describing 
Green's Function Monte Carlo simulations of light nuclei and neutron matter based on AVI 8 



32 



33 



34], 



35 



37 



38, 



as well as other phenomenological potentials 27, 

rt 

review article detailing this method [4j as well as a more recent set of lecture notes |40j 



391 ] . There is a 
A 



related technique known as auxiliary-field diffusion Monte Carlo simplifies the spin structure 



of the same calculations by introducing auxiliary fields 



41 



42 



43 



44 



45 



46|. The No-Core 



Shell Model (NCSM) is a different approach to light nuclei which produces approximate 
eigenvectors in a reduced vector space. There have been several NCSM calculations using 
various different phenomenological potential models [47J, |48|, There are also NCSM 
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calculations which have used nuclear forces derived from chiral effective field theory [5ll. 152. 
531 ] . Recently there has also been work in constructing a low-energy effective theory within 
the framework of truncated basis states used in the NCSM formalism 54( . A benchmark 



comparison of many of the methods listed above as well as other techniques can be found 
in Ref. [55]. 



In this article we describe recent work by several different collaborations which combine 
the framework of effective field theory with computational lattice methods. The idea of 
lattice simulations using effective field theory is rather new. The first quantum lattice study 



of nuclear matter appears to be Ref. 



56] , which used a momentum lattice and the quantum 



hadrodynamics model of Walecka 57] . The first study combining lattice methods with an 



effective theory for low-energy nuclear physics was Ref. [58(. This study looked at infinite 
nuclear and neutron matter at nonzero density and temperature. After this there appeared 
a computational study of the attractive Hubbard model in three dimensions {59I, as well as 
a paper noting the absence of sign oscillations for nonzero chemical potential and external 
pairing field 60]. Another study looked at nonlinear realizations of chiral symmetry with 



static nucleons on the lattice 



611 ] . and there were also a number of investigations of chiral 



perturbation theory with lattice regulation SIHBi. Th, was Mowed by the Sr S t 
many-body lattice calculation using chiral effective field theory [65J. From about this time 
forward there were a number of lattice calculations for cold atoms and low-energy nuclear 
physics which we discuss in this article. 

The lattice effective field theory approach has some qualitative parallels with digital 
media. In digital media input signals are compressed into standard digital output that can 
be read by different devices. In our case the input is low-energy scattering data, and the 
digital format is effective field theory defined with lattice regularization. The process of 
sampling and compression consists of matching low-energy scattering data using effective 
interactions up to some chosen order in power counting. By increasing the order, the 
accuracy in describing low-energy phenomena can be systematically improved. 

Just as standard digital format enables communication between different devices, lattice 
effective field theory enables the study of many different phenomena using the same lattice 
action. This includes few- and many-body systems as well as ground state properties and 
thermodynamics at nonzero temperature and density. Another attractive feature of lattice 
effective field theory is the direct link with analytic calculations using effective field theory. 
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It is straightforward to derive lattice Feynman rules and calculate diagrams using the same 
theory used in non-perturbative simulations. At fixed lattice spacing all of the systematic 
error is introduced up front when defining the low-energy lattice effective field theory and 
not determined by the particular computational scheme used to calculate observables. This 
allows for a wide degree of computational flexibility, and one can use a number of efficient 
lattice methods already developed for lattice QCD and condensed matter applications. This 
includes cluster algorithms, auxiliary-field transformations, pseudofermion methods, and 
non-local configuration updating schemes. We discuss all of these techniques in this article. 
We also review the relevant principles of effective field theory as well as different formalisms 
and algorithms used in lattice calculations. Towards the end we discuss some recent results 
and compare with results obtained using other methods. 



II. EFFECTIVE FIELD THEORY 



Effective field theory provides a systematic approach to studying low-energy phenomena 
in few- and many-body systems. We give a brief overview of the effective range expansion 
and the application of effective field theory to cold atoms and low-energy nuclear physics. 
A more thorough review of effective field theory methods applied to systems at nonzero 



density can be found in Ref. 



66|. 



A. Effective range expansion 

At sufficiently low momentum the cross-section for two-body scattering is dominated by 
the S-wave amplitude, and higher partial waves are suppressed by powers of the relative 
momentum. The S-wave scattering amplitude for two particles with mass m and relative 
momentum p is 

MP) = TT (2-1) 

m p cot do — ip 

where 5q is the S'-wave phase shift. At low momentum the S'-wave phase shift for two- 
body scattering with short-range interactions can be written in terms of the effective range 



expansion 



67|, 



p cot 5 = — + \r cS p 2 + ■■■ . (2.2) 

^scatt ^ 
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Here a sca tt is the S'-wave scattering length, and r e g is the S'-wave effective range. The radius 
of convergence of the effective range expansion is controlled by the characteristic length scale 
of the interaction. For example in low-energy nuclear physics the range of the two-nucleon 
interaction is set by the Compton wavelength of the pion. The generalization of the effective 
range expansion to partial wave L has the form 

p 2L+l cot 5 L = -— + \r L p 2 + ■■■ . (2.3) 
a L 2 

The 5l phase shift scales as 0{p 2L+l aL) in the low-momentum limit, and higher-order terms 
are suppressed by further powers of p 2 . This establishes a hierarchy of low-energy two-body 
scattering parameters for short-range interactions. For particles with intrinsic spin there is 
also some mixing between partial waves carrying the same total angular momentum. 

For many interacting systems we can characterize the low-energy phenomenology ac- 
cording to exact and approximate symmetries and low-order interactions according to some 
hierarchy of power counting. This universality is due to a wide disparity between the 
long-distance scale of low-energy phenomena and the short-distance scale of the underlying 
interaction. In some cases the simple power counting of the effective range expansion must 
be rearranged or resummed in order to accommodate non-perturbative effects. We discuss 
this later in connection with singular potentials and three-body forces. A recent review of 



68] 



universality in few-body systems at large scattering length can be found in Ref. 

In many-body systems a prime example of universality is the unitarity limit. The 
unitarity limit describes attractive two-component fermions in an idealized limit where the 
range of the interaction is zero and the scattering length is infinite. The name refers to the 
fact that the .S'-wave cross-section saturates the limit imposed by unitarity, <Jq{p) — 47r/p 2 , 
for low momenta p. While the unitarity limit has a well-defined continuum limit and strong 
interactions, at zero temperature it has no intrinsic physical scale other than the interparticle 
spacing. 

Phenomenological interest in the unitarity limit extends across several subfields of physics. 
The ground state of the unitarity limit is known to be a superfluid with properties in between 
a Bardeen-Cooper-Schrieffer (BCS) fermionic superfluid at weak attractive coupling and a 



Bose-Einstein condensate (BEC) of bound dimers at strong attractive coupling 69|, |7CJ, 



711 ] . It has been suggested that the crossover from fermionic to bosonic superfluid could 



be qualitatively similar to pseudogap behavior in high-temperature superconductors 



72] 



In nuclear physics the unitarity limit is relevant to the properties of cold dilute neutron 
matter. The neutron scattering length is about —18.5 fm while the range of the interaction is 
comparable to the Compton wavelength of the pion, m" 1 « 1.4 fm. Therefore the unitarity 
limit is approximately realized when the interparticle spacing is about 5 fm. Superfluid 



neutrons at around this density may exist in the inner crust of neutron stars 73j, 1 741 ]. 



B. Effective field theory for cold atoms 



Physics near the unitarity limit has been experimentally observed in cold degenerate 
gases of 6 Li and 40 K atoms. Alkali atoms are convenient for evaporative cooling due to 
their predominantly elastic collisions. For sufficiently dilute gases the effective range and 
higher partial wave effects are negligible whi 
a magnetically-tuned Feshbach resonance 



Feshbach resonances can be found in Ref. 179. 




e the scattering length can be adjusted using 
Overviews of experiments using 



801 ] , and there are a number of reviews covering 



81 



82]. 



the theory of BCS-BEC crossover in cold atomic systems 

At long distances the interactions between alkali atoms are dominated by the van der 
Waals — Ce/r 6 interaction. Power-law interactions complicate the effective range expansion 
by producing a branch cut in each partial wave at p = 0. For the van der Waals interaction 
the expansion in p 2 is an asymptotic expansion coinciding with the effective range expansions 



in Eq. ( 12. 2ft and (12. 3ft up through terms involving a scat t, r e g, and a\ |83|, |84j. Beyond this the 
asymptotic expansion involves powers of p 2 times hap 2 or odd powers of p. All of the work 
discussed in this article involves low-energy phenomena where these non-analytic terms can 
be neglected. 

The low-energy effective field theory for the unitarity limit can be derived from any 
theory of two-component fermions with infinite scattering length and negligible higher-order 
scattering effects at the relevant low-momentum scale. For example the two fermion compo- 
nents may correspond with dressed hyperfine states \ f,m,f) = 1 9/2, — 9/2) and 1 9/2, — 7/2) 
of 40 K with interactions given either by a full multi-channel Hamiltonian or a simplified 



two-channel model 



85 



1861 . |87j . The starting point does not matter so long as the S- wave 
scattering length is tuned to infinity to produce a zero-energy resonance. 

In our notation m is the atomic mass and dj and a\ are annihilation and creation operators 
for two hyperfine states. We label these as up and down spins, i =|, j, even though the 
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connection with actual intrinsic spin is not necessary. We enclose operator products with 
the symbols :: to indicate normal ordering, where creation operators are on the left and 
annihilation operators are on the right. The effective Hamiltonian at leading order (LO) is 

HhO = Hfree + Vlo, (2.4) 

where 




(2.6) 

and p aT > a (f) is the particle density operator, 

/' a (r1 = £ a\{r)a t {r). (2.7) 

The coefficient C depends on the cutoff scheme used to regulate ultraviolet divergences 
in the effective theory. Higher-order effects may be introduced systematically as higher- 
dimensional local operators with more derivatives and/or more local fields. 

C. Pionless effective field theory 

For nucleons at momenta much smaller than the pion mass, all interactions produced by 
the strong nuclear force can be treated as local interactions among nucleons. The effective 
Hamiltonian in Eq. (12.41) also describes the interactions of low-energy neutrons at leading 
order. For systems with both protons and neutrons we label the nucleon annihilation 
operators with two subscripts, 



&o,o — cl-\, p , ao,i — Q>i,m (2-8) 

a i,o = a>i,p, a i,i = cti,n- (2-9) 

The first subscript is for spin f, j and the second subscript is for isospin p, n. We use as 
with S = 1 , 2, 3 to represent Pauli matrices acting in spin space and 17 with I = 1 , 2, 3 to 



represent Pauli matrices acting in isospin space. The same letters S and I are also used 
to indicate total spin and total isospin quantum numbers, but the intended meaning will 
be clear from the context. If we neglect isospin breaking and electromagnetic effects, the 
effective theory has exact SU(2) spin and SU(2) isospin symmetries. 

9 



LO 



d s r 



P a ' a (f) 



Let us define the total nucleoli density 



P aV (r*) = 



(2.10) 



The total nucleon density is invariant under Wigner's SU(4) symmetry mixing all spin and 



isospin degrees of freedom 



881 ] . Using as and rj, we also define the local spin density, 



isospin density, 



and spin-isospin density, 



Ps ( f ) = a h ( f ) fa] «' ai ' >i ^ ' 

i,j,i'=0,l 



mj'=o,i 



P a s?(r) = Yl a h^ ^ s ^i' ( T/ W 

i,j,i',j'=0,l 



a i',j'{ r ) 



At leading order the effective Hamiltonian can be written as 



(2.11) 



(2.12) 



(2.13) 



Hlo — H ircc + VlO) 



(2.14) 



where 



free 



— V 

2m ^— ' 



dVVa^.(f) ■ Vaij(f), 

»,i=o,i'' 

^lo = V + V P + V {3N \ 

2 



V — ^ / d''r : 



7=1,2,3 

D 
"6 



1 2 



y (37V) = ^ / d 3 f 



P at ' a (r)l 3 



Due to an instability in the limit of zero-range interactions 89j, the SU(4)-symmetric three 



(2.15) 

(2.16) 
(2.17) 

(2.18) 
(2.19) 



90 



91 



92] 



nucleon force V^ 3N ^ is needed for consistent renormalization at leading order 
With the constraint of antisymmetry there are two independent S'-wave nucleon-nucleon 
scattering channels. These correspond with spin-isospin quantum numbers S = 1, / = 

and S — 0, I — I. i Some analytic methods used in pionless effective field theory are 

discussed in Ref. 



93 



94j. A general overview of methods in pionless effective field theory 



can be found in recent reviews 



95 



96 



97] - 



10 



D. Chiral effective field theory 



For nucleon momenta comparable to the pion mass, the contribution from pion modes 
must be included in the effective theory. In the following q denotes the t-channel momentum 
transfer for nucleon-nucleon scattering while k is the w-channel exch ang ed momentum trans- 



98 



99j the nucleon-nucleon 



fer. At leading order in the Weinberg power-counting scheme 
effective potential is 

Hlo = H {rcc + Vio, (2.20) 

V LO = V + Vl2 + V orvr. (2.21) 

# free , V, Vp are defined in the same manner as in Eq. fT2~T5l) . ([2TTT|) . f[2TT8]) . V OPEP is the 
instantaneous one-pion exchange potential, 

^° PEP = E / d ^id^G Sl s 2 (ri - r 2 ) : p&Jft) P&Jft) (2-22) 

Si,S 2 ,I=l,2,3^ 

where the spin-isospin density p^f is defined in Eq. (12 . 13[) and 

For our physical constants we take m = 938.92 MeV as the nucleon mass, m n = 138.08 MeV 
as the pion mass, f n = 93 MeV as the pion decay constant, and Qa = 1-26 as the nucleon 
axial charge. 

The terms in Vlo can be written more compactly in terms of their matrix elements with 
two-nucleon momentum states. The tree-level amplitude for two-nucleon scattering consists 
of contributions from direct and exchange diagrams. However for bookkeeping purposes we 
label the amplitude as though the two interacting nucleons are distinguishable. We label 
one nucleon as type A, the other nucleon as type B, and the interactions include densities 
for both A and B. For example the total nucleon density becomes 

p a\a _^ p a\,a A + p a< B ,a B _ (2.24) 



The amplitudes are then 



A(V) = C, (2.25) 
A(V l2 )=CpJ2rfrf, (2.26) 
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1 j " V2/J ? 2 + < ' { } 

At next-to-leading order (NLO) the effective potential introduces corrections to the two 

LO contact terms, seven independent cont act terms carrying two p owers of momentum, and 



instantaneous two-pion exchange (TPEP) [100 



101 



102 



103 



1041 ]. We write this as 



^nlo = V LO + AV<°> + W + ^™\ 



The tree-level amplitudes for the new contact interactions are 



A(AV) = AC, 



i 

A(V q2 )=C q2 q 2 , 

A(V l2ig 2)=C l2 , q2 q 2 J2 T ? T ? 

I 

A {V S 2, q 2) = C S 2^q 2 

s 

A (Ks 2 ,/ 2 ,? 2 ) = CW 2 ,<7 2 <? 2 a s a s t i t i > 

S I 
S S' 

A (Vp,( ff .s)a) = C P>{q . s) 2 Y T i T f qsa s ^ ,cr f'' 

I S S' 

A (V {iqxS yk) = iC( iqxS ).k ^ eisi'Qi {& A + <? B ) s h>- 

i,s,v 



The amplitude for NLO two-pion exchange potential is 



105 



106] 



^(C P ) 



384vr 2 /4 
3^ 4 



Ami {$9 a ~ 49a ~l)+q 2 (23<?1 - 10<& - l) 



where 



2g \/4m 2 + g 2 - g 



Recent reviews of chiral effective field theory can be found in Ref. 



95 



96 



97]- 



(2.28) 



(2.29) 

(2.30) 

(2.31) 
(2.32) 

(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 

4:8g\ml ' 
4m 2 + q 2 

(2.38) 
(2.39) 
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E. Three-nucleon forces 



The systematic framework provided by effective field theory becomes very useful when 
discussing the form of the dominant three-nucleon interactions. Few-nucleon forces in chiral 
effect ive f ield theory beyond two nucleons were first discussed qualitatively in Ref. 99]. In 



Ref. 



1071 ] it was shown that the three-nucleon terms at NLO cancelled, and the leading 



three-nucleon effects appeared at next-to-next-to leading order (NNLO) in Weinberg power 
counting. 

The NNLO three-nucleon effective potential arises from a pure contact potential, V^, 



contact 5 



one-pion exchange potential, Vqpe , and a two-pion exchange potential, . Parts of the 



(3N) 



NNLO three-nucle on potentia 
nucleon potentials 



26 



27. 



108 



are also cont ained in a number of phenomenological three- 



109 



110 



llll ] . However there is clear value in identifying the 



full set of leading interactions. Similar to our description above for two-nucleon scattering, 
we write the tree-level amplitude for three-nucleon scattering where the first n ucleon is of 
type A, the second nucleon type B, and the three type C. The amplitudes are [112], 1 1 1 3 ] 



A V r 



r(3N) \ 

' contact J 



D. 



contact i 



AY 



.(37V) 



OPE 



Jn perm A,B,C 



q\ + ml 



^2.40) 



(2.41) 



A I V (3N) 

./i I 'TPE 



C 3 



9l 



E 



(qa ■ &a) (qb ■ o B ) {qa ■ Qb) 



- c r 



Ta ■ &a) (qb ■ °b) 



^ pe J^, B ,c(ii + <niB + mi) 



[TA ■ TB) 



9a ST (U ■ a a) {Qb ■ gg) 

° 4 8 ft ^ (ql + ml) (q% + ml) 

J7T perm A,B,C KHA kIWB ttJ 



[(q A x q B ) ■ a c ] [{t a x t b ) • t c ] ■ (2.42) 



In our notation q B , qc are the differences between final and initial momenta for the 
respective nucleons. The summations are over permutations of the bookkeeping labels 
A, B, C. 

The coefficients 01,3,4 are irirNN interaction t erm s in the chiral Lagrangian and are de- 
termined from fits to low-energy scattering data J 
-Dcontact and -Dope are cutoff dependent. In Ref. 



4L The remaining unknown coefficients 



1131 ] these were fit to the triton binding 



energy and spin-doublet neutron-deuteron scattering length. The resulting NNLO effective 
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potential was shown to give a prediction for the isospin-symmetric alpha binding energy 
accurate to within a fraction of 1 MeV. 



F. Non-perturbative physics and power counting 

When non-perturbative processes are involved, reaching the continuum limit and power 
counting in effective field theory can sometimes become complicated. The two-component 
effective Hamiltonian for cold atoms introduced in Eq. ( 12.41) has no such complications. 
Ultraviolet divergences can be absorbed by renormalizing the interaction coefficient C, and 
the cutoff momentum can be taken to infinity. Similarly the leading-order pionless effective 
Hamiltonian in Eq. (I2.14p has a well-defined continuum limit if we neglect deeply-bound 
three-body states that decouple from the low-energy effective theory. While these deeply- 
bound states generate instabilities in numerical simulations they can be removed by hand 



in semi-analytic calculations 90l. l9lL 192]. 



In chiral effective field theory there has been considerable study on the consistency of 
the Weinberg power counting scheme at high momentum cutoff. Complications arise from 
the singular behavior of the one-pion exchange potential. In order to avoid unsubtracted 



ultraviolet divergences produced by infinite iteration of the one-pion exchange potentia 



alte rnative scheme was proposed where pion exchange is treated perturbatively 



115 



an 



U6, 



1171 ] . This approach, KSW power counting, allows for systematic control of the ultraviolet 



divergence structure of the effective theory. Unfortunately the convergen ce a t higher order 



118] 



is poor in some partial waves for momenta comparable to the pion mass 

The most divergent short-distance part of the one-pion exchange potential is a f~ 2 r~ 3 sin- 
gularity arising from the tensor force in the spin-triplet channel. There are also subleading 
divergences at r = which contain explicit factors of the p ion m ass. Based on this observa- 



tion another power counting scheme was proposed in Ref. [119]. This new scheme coincides 
with KSW power counting in the spin-singlet channel. But in the spin-triplet channel the 
most singular piece of the one-pion exchange potential is iterated non-perturbatively, while 
the rest is incorporated as a perturbative expansion around = 0. 



More recently a different power counting modification was proposed in Ref. 1201 ] . In 
this approach the one-pion exchange potential is treated non-perturbatively in lower angular 
momentum channels along with higher-derivative counterterms promoted to leading order. 
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These counterterms are used to cancel cutoff dependence in channels where the tensor force is 
attractive and strong enough to overcome the centrifugal barrier. Advantages over Weinberg 
power counting at leading order were shown for cutoff momenta much greater than the pion 



mass. Further investigations of this approach i n hig ! 
with one-pion exchange were considered in Ref. 121 



rer partial waves and power counting 



122 1 . 



The choice of cutoff momentum and power counting scheme in lattice effective field theory 
is shaped to a large extent by computational constraints. For two-nucleon scattering in chiral 
effective field theory, small lattice spacings corresponding with cutoff momenta many times 
greater than the pion mass are no problem. However at small lattice spacing significant 
numerical problems appear in simulations of few- and many-nucleon systems. In attractive 
channels one must contend with spurious deeply-bound states that spoil Euclidean time 
projection methods (a technique described later in this review). In channels where the 
short-range interactions are repulsive a different problem arises. In auxiliary-field and 
diagrammatic Monte Carlo (methods we discuss later in this review), repulsive interactions 
produce sign or complex phase oscillations that render the method ineffective. Due to 
these practical computational issues one must settle for lattice simulations where the cutoff 
momentum is only a few times the pion mass, and t he ad vantages of the improved scheme 
over Weinberg power counting are numerically small 123 ]. 



III. LATTICE FORMULATIONS FOR ZERO-RANGE ATTRACTIVE TWO- 
COMPONENT FERMIONS 

In this section we introduce a number of different lattice formulations using the example 
of zero-range attractive two-component fermions described by i?LO i n Eq. ( 12. 4ft . In Fig. (CQ) 
we show a schematic diagram of the different lattice formulations. The numbered arrows 
indicate the discussion order in the text. 

Throughout our discussion of the lattice formalism we use dimensionless parameters and 
operators corresponding with physical values multiplied by the appropriate power of the 
spatial lattice spacing a. In our notation the three-component integer vector n labels the 
lattice sites of a three-dimensional periodic lattice with dimensions L 3 . The spatial lattice 
unit vectors are denoted I — 1, 2, 3. We use n t to label lattice steps in the temporal 
direction, and L t denotes the total number of lattice time steps. The temporal lattice 
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FIG. 1: A schematic diagram of different lattice formulations. The numbered arrows indicate the 
discussion order in the text. 



spacing is given by a t , and a t = a t /a is the ratio of the temporal to spatial lattice spacing. 
We also define h = a t /(2m), where m is the fermion mass in lattice units. 

A. Grassmann path integral without auxiliary field 

For two-component fermions with zero-range attractive interactions we start with the lat- 
tice Grassmann path integral action without auxiliary fields. It is the simplest formulation 
in which to derive the lattice Feynman rules. Hence it is usef ul for bo th analytic lattice 



calculations and diagrammatic lattice Monte Carlo simulations [124J . Il25 |. 

We let Cj and c* be anticommuting Grassmann fields for spin i =|, [. The Grassmann 
fields are periodic with respect to the spatial lengths of the L 3 lattice, 

Cj(n + LI, n t ) = Ci(n + L2, n t ) = Cj(n + L3, n t ) = c^n, n t ), (3.1) 

c*(n + Ll } n t ) = c*(n + L2,n t ) = c*(n + L3,n t ) = c*(n,n t ), (3.2) 
and antiperiodic along the temporal direction, 

Cj(n, rit + L t ) = -Cj(n, n t ). (3.3) 

c*(n,n t + L t ) = -c*(n,n t ). (3.4) 
We write DcDc* as shorthand for the integral measure, 

DcDc* = | [ dci(n,n t )dc*(n,n t ). (3.5) 

n,n t ,i=1,l 
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We use the standard convention for Grassmann integration, 

dci(n,n t ) = / dc*(n,n t ) = 0, (3.6) 



dci(n,nt)ci(n,n t ) = J dc*(n,n t )c*(n,n t ) = 1 (no sum on i). (3.7) 

Local Grassmann densities p-\, Pi, p are defined in terms of bilinear products of the Grass- 
mann fields, 

p^{n,n t ) = c*(n,n t )c^(n,n t ), (3.8) 
Pl(n,n t ) = clfan^Cifant), (3.9) 
p(n, n t ) = p T (n, n t ) + p^n, n t ). (3.10) 
We consider the Grassmann path integral 

Z = J DcDc*exp[-S (c,c*)], (3.11) 

where 

S(c,c*) = S iree (c,c*) + Cat^2pi(n,nt)pi(n,n t ). (3.12) 

n,n t 

The action S(c, c*) consists of the free nonrelativistic fermion action 

Sfree(c,c*)= [c*(n, n t )ci(n, n t + 1) - (1 - 6/1) c£(n, n t )cj(n, n f )] 

- ^ jc*(n,n t )c;(n + Z,n t ) + c*(n,nt)ci(n - l,n t ) , (3.13) 

n,n t ,i=T,J, (=1,2,3 

and a contact interaction between up and down spins. We consider the case where the coef- 
ficient C is negative, corresponding with an attractive interaction. Since we are considering 
nonrelativistic lattice fermions with a quadratic dispersion relation, the lattice doubling 
problem associated with relativistic fermions does not occur. 

In the grand canonical ensemble a common chemical potential p is added for all spins. 
In this case the /i-dependent path integral is 

Z(fi) = J DcDc*exp{-S(c,c*,p)}, (3.14) 

where 

S(c,c*,p) = S(e^c,c*) + Yl l(l-e^)c*(n,n t )c t (n,n t + l)], (3.15) 

n,n t ,i=1,l 

and S(e fiat c, c*) is the same as S(c, c*) defined in Eq. (13.121) . but with c replaced by e mt c. 
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B. Transfer matrix operator without auxiliary field 



Let a and denote fermion annihilation and creation operators satisfying the usual 
anticommutation relations 

{a, a} = {a\a ] } = 0, (3.16) 

{a,a i } = l. (3.17) 



For any function / {a) , a) we note the identity 126 1 



Tr[: f (a\a) : ] = J dcdc* e 2c * c f (c* , c) , 



(3.18) 



where c and c* are Grassmann variables. As before the :: symbols in Eq. (I3.18P indicate 
normal ordering, and the trace is evaluated over all possible fermion states. This result can 



be checked explicitly using the complete set of possible functions {l 



, a, a\ a) a 



}■ 



It is useful to write Eq. (I3.18P in a form that resembles a path integral over a short time 
interval with antiperiodic boundary conditions, 



Tr [: / (at, a) :]= J c/c(0)rfc*(0)e c * (0)[c(0) - c(1)1 / [c*(0), c(0)] , 



c(l) = -c(0). 



(3.19) 
(3.20) 



This result can be generalized to products of normal-ordered functions of several creation 
and annihilation operators. Let a^(n) and a\(n) denote fermion annihilation and creation 
operators for spin i at lattice site n. We can write any Grassmann path inte gral with 
instantaneous interactions as the trace of a product of operators using the identity [1261 . Il27 ] 



Tr 







'Lt-l 



a\{n'),ai{n)\ : x • • • x : F a],{n'),ai{n) 



j DcDc* exp < ^ ^ c*(n, n*) [c;(n, n*) - Cj(n, n t + 1)] > 

I nt=0 n,i J 



Lt-l 



X 



II ^n t [4(^,nt),Cj(n,nt)] 



(3.21) 



n t =0 



where Ci(n,L t ) = — q,(n, 0). 

Let us define the free nonrelativistic lattice Hamiltonian 



free 



m 



£ 4( 



n aAn - 



n,i=t,| 



1 

2m 



[aI(n)oi(n + l) + at(n)a;(n - f)l , (3.22) 



n,i=T,4. Z=l,2,3 
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as well as the lattice density operators 

pf a (n) = a{(n)a T (n), (3.23) 

pf a (n) = a\(n) ai (n), (3.24) 

p ata (n) = pf a {n) + pf a {n). (3.25) 

Using the correspondence Eq. (13.211) . we can rewrite the path integral Z defined in Eq. (13. lip 
transfer matrix partition function, 

Z = Tr (M Lt ) , (3.26) 

where M is the normal-ordered transfer matrix operator 



M =: exp 



H {ree a t -Ca t y2 P f a (n)pf a (n) 



(3.27) 



Roughly speaking the transfer matrix operator is the exponential of the Hamiltonian oper- 
ator over one Euclidean lattice time step, e~ Hat . In order to satisfy the identity Eq. (13.211) . 
we work with normal-ordered transfer matrix operators. In the limit of zero temporal lattice 
spacing, a t — > 0, we obtain the Hamiltonian lattice formulation with Hamiltonian 

H = tf free + CY,pf a (n)pf a (n). (3.28) 

n 

This is also the defining Hamiltonian for the attractive Hubbard model in three dimensions. 

In the grand canonical ensemble the effect of the chemical potential is equivalent to 
replacing M by 

M(ji) = M exp //id* P aU ^ | • ( 3 - 29 ) 

For the Hamiltonian lattice formulation the effect of the chemical potential has the familiar 
form 

H(p) = H iTee + Cj2pf a (n)Pi a (n) ~ j"X> ata ^)- ( 3 - 3 °) 

it n 

C. Grassmann path integral with auxiliary field 

We can re-express the Grassmann path integral using an auxiliary field coupled to the 
particle density. This lattice formulation has been used in several lattice studies at nonzero 
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temperature [60|, |65|, 1 12a . Il29l . Il30l . Il3ll . 1 1321 . 1 1 33f ] . Due to the simple contact interaction 
Pl(n,nt)pi(n,n t ) and the anticommutation of Grassmann variables, there is a large class of 
auxiliary-field transformations which reproduce the same action. 

Let us write the Grassmann path integral using the auxiliary field s, 



where 



Z = Y\ J d A s(n,n t ) J DcDc* exp [— S A (c, c*, s)\ 

n,nt 



S A (c, c*, s) = S bee (c, c*) - ^ A [ s (™> n *)] n *) 



(3.31) 



(3.32) 



n,n t 



One possible example is a Ga ussian- int egral transformation similar to the original Hubbard- 



Stratonovich transformation 



134 



d A s(n, n t 



135| where 

i r+oc 



ds(n, n t )e~ 



\s 2 (n,n t ) 



V27T j -oo 

A[s(n,n t )] = \J-Cot t s{fi, n t ). 



(3.33) 



(3.34) 



Anot 
Ref. 



rer possibility is a discrete auxiliary- field transformation similar to that used in 



1361 ] . In our notation this can be written as 

1/2 



f d A s(n,n t ) = I c 

J J -1/2 



ds(n, n t ) 



A [s(n, nt)] = \f-Ca t sgn [s( 



n,n t 



(3.35) 
(3.36) 



1371 ] the 



where sgn equals +1 for positive values and —1 for negative values. In Ref. 
performance of four different auxiliary-field transformations were compared. 

We intentionally leave the forms for d A s(n,nt) and A [s(n,rit)] unspecified, except for a 
number of conditions needed to recover Eq. (13.111) upon integrating out the auxiliary field 
s. The first two conditions we set are 



d A s{n, n t )l 



d A s(n, n t ) A [s(n, n t )] = 0. 



(3.37) 



(3.38) 



Since all even products of Grassmann variables commute, we can factor out the term in 
Eq. (13.311) involving the auxiliary field s at n,nt- To shorten the notation we temporarily 
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omit writing n, n t explicitly. We find 

J d A s exp [ A (s) (p T +Pi)] = J d A s [1+ A (s) (p T + Pl ) + A 2 (s) P]Pl ] 

= 1+ / d A s A 2 (s) p T p| = exp J d A s A 2 (s) P i Pl 



Therefore the last condition needed to recover Eq. ( 13. lip is 

- Ca t = / d A s A 2 (s) . 



(3.39) 



(3.40) 



In the grand canonical ensemble, the auxiliary-field path integral at chemical potential p is 

= II J d A s(n,n t ) J DcDc* exp[-S A (c,c*,s, p)] , (3.41) 



where 



S A (c, c*, s, p) = S A (e^c, c*,s)+ Yl K 1 - e ^ c ^ n *) c ^' n * + X )] • ( 3 - 42 ) 



n,nt,i=],l 



D. Transfer matrix operator with auxiliary field 



Using Eq. (I3.2ip and (13.311) we can write Z as a product of transfer matrix operators 
which depend on the auxiliary field, 



n,n t 



d A s(n, n t ) 



Tr{M A (s,L t -l) M A (s,0)} 



where 



This form has 



M A (s, m) = : exp j-# frce a 4 + J2 A ^ nt ^ ^"( ^j 



(3.43) 



(3.44) 



3een used in a number of lattice simulations 



141 



142 



143 



58 



132 



133 



137 



138 



139 



144 



145 



146 



HO, 



1471 ] . In some of these studies the Hamiltonian limit a t — > 



is also taken. 

In the grand canonical ensemble at chemical potential p the partition function is 



Z( P ) = U\j d 

n,n t 



l A s{n,n t 



Tr {M A (s, L t - 1, p) M A (s, 0, p)} , (3.45) 



where M A (s,n t , p) is defined as 



M A (s, n t , p) = M A (s, nt) exp 



|pa^p ata (n)| 



(3.46) 
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E. Improved lattice dispersion relations 



In Ref. 138l . 1 1461 . 1 1471 ] the transfer matrix operator at chemical potential \i was written 



as 



exp 



~ ( #lr.v - /'-V 



exp 



n 



exp 



~ ( //iv,,. - /'.V 



(3.47) 



in the Hamiltonian limit, at — > 0, with 

N = Y,P ata (n)- 



(3.48) 



This is different from M(/z) in Eq. (I3.29p . but the two are the same in the Hamiltonian limit. 
The exponential interaction term in Eq. (13.471) was treated using a discrete auxiliary field. 
Also the matrix elements of Hf ree were computed by Fast Fourier Transform in momentum 
space using the quadratic dispersion relation 



CO 



(quad) 



(p) 



- E 

2m ^ 



Pi: 



(3.49) 



«=1,2,3 



with pi defined in the first Brillouin zone, \pi\ < it. The motivation for this approach was 
to remove errors associated with the standard lattice dispersion relation 



m A — ' 



cos pi) 



(3.50) 



J=l,2,3 



In Ref. 



130 



13l| lattice calculations at nonzero temperature and large scat terin g length 



found significant errors due to lattice artifacts. A detailed analysis in Ref. [14a ] showed 
that the large errors were produced by broken Galilean invariance on the lattice. As an 
alternative to the momentum space approach in Eq. (13.491) . improved lattice dispersions 
were investigated that could be derived from local lattice actions. 

A class of improved single-particle dispersion relations can be defined on the lattice, 

uj(n) ^ = ^ E E (-i)S- n)cos fe)- (3.5i) 

3=0,1,2,- 1=1,2,3 

uj(°\p) corresponds with the standard action, uj^ip) is the 0(a 2 )-improved action, and so 
on. The improved actions eliminate lattice artifacts in the Taylor expansion of a; < - n ^(p) about 
p = 0, 



2n+2 



(3.52) 



2=1,2,3 
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The lattice action corresponding with u^ n ' contains hopping terms in each spatial direction 

(ri) 

that extend n lattice steps beyond the nearest neighbor. The hopping coefficients for 
actions up to 0(a 4 ) are shown in Table [fl 

TABLE I: Hopping coefficients for lattice actions up to 0(a 4 ). 





standard 


0(a 2 )-improved 


0(a 4 )-improved 


^0 


1 


5 
4 


49 

36 


Vl 


1 


4 
3 


3 
2 


V2 





1 

12 


3 

20 


V3 








1 

90 



In addition to these improved actions, new lattice actions called well-tempered actions 
were also introduced. These were defined implicitly in terms of their dispersion relation, 

(3.53) 



J wtn \p) = u {n ~ l \p) + c [u {n \p) - uj {n - l) ( 
where the unknown constant c was determined by the integral constraint, 



7T 7T 7T 




dp 1 dp 2 dp 3 



-7T — IX — 7T 



(wtn) . 



2=1,2,3 



0. 



(3.54) 



At nonzero temperature and large scattering length the local well-tempered action corre- 
sponding with cj( wtl ) was shown to be comparable in accuracy to the nonlocal action defined 
by w (quad) [l48[. 



IV. LATTICE FORMULATIONS FOR LOW-ENERGY NUCLEONS 



A. Pionless effective field theory 



Analogous with the continuum densities in Eq. (12.101) . (12. lip . (12.121) . and (12 . 13j) . we define 
the lattice operators 

" ' (ft), (4.1) 



i,j=0,l 



Ps 



n 



alA^^sli'diijin), 



(4.2) 



ij',»'=0,l 
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■An) 



pf' a (n)= Nii'°iJ 

ij',j'=0,l 

Psf(n)= J2 a h^ ^sluf [ r i)jf a i',j'(n). 

i,j,i',j'=Q,l 

At leading order in pionless effective field theory, 



(4.3) 
(4.4) 



where 



M 



Z = Tr (M Lt ) , 



■■: expl -H {ree a t - -Ca t ^ \p a ^' a {n) 



(4.5) 



ft, I n 



n 



(4.6) 



This formalism was used to study the triton and three-body forces on the lattice 149| . The 
triton can be regarded as an approximate example of the Efimov effect, which in the limit of 
zero range and infinite s cattering length predicts a geometric sequence of trimer bound states 



68. 



90 



91 



92 



150 



15 1| . The Efimov effect is not possible for two-component fermions due 



to Pauli exclusion but is allowed for more than two components. Once the binding energy 
of the trimer s ystem is fixed, the binding energy of the four-body system is also determined 



152 



153 



1541 ] . This is in analogy with the Tjon line relating the nuclear binding energies 



of 3 H and 4 He. In two dimensions a different geometric sequence has been predicted for 
zero-range attractive interactions. In this case the geometric sequen ce describes th e binding 



energy of iV-body clusters as a function of N in the large N limit 



155 



156 



1571 ] . These 



two-dimensional clusters have been studied using lattic e eff ective field theory for up to 10 



particles and the geometric scaling has been confirmed 



i3- 



B. Pionless effective field theory with auxiliary fields 



In terms of auxiliary fields 



n,n t 



d A s(n, n t ) 



n 

n,nt,I 



xTr{M A (s,L t -l) 



'27T J-oo 

■M a (s,0)} 



(4.7) 
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where the auxiliary-field transfer matrix is 



M A (s, n t ) = : exp I -H {ree a t + ^ A [s(n, n t )} p ata (n) 

I n 

+ i^/cja t ^2s 1 (n,n t )pf a { 



(4.8) 



Let (^4 fc ) be the expectation value of the k th power of A with respect to the measure o^s, 

<A fe ) = Jd A s(n,n t ) {A[s(n,n t )]} k , k — 0, 1,2,3,4. (4.9) 
In order to reproduce the interactions in Eq. ( 14. 6 j) we require that 



(A°> = 1, <A 1 > = 0, <A 2 ) = -Ca t , <A 3 > 



(A 4 ) = 3C 2 a 2 . 

(4.10) 

The existence of a positive definite measure g^s and real- valued A is essential for Monte 
Carlo simulations without sign and phase oscillations. Sufficient and necessary conditions 
for the existence of a positive definite <1as and real-valued A is known in the mathematics 
iterature as t he truncated Hamburger moment problem. This problem has been solved 
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1601 ] . and the conditions are satisfied if and only if the block-Hankel matrix, 



(A°) {A 1 } (A 2 ) 
(A 1 ) (A 2 ) (A 3 ) 
(A 2 ) (A 3 ) (A 4 ) 



1 





-Ca t 
-Ca t —Da t 



-Ca t -Dot 3C 2 a 2 t 



D 2 al 



(4.11) 



With an 



is positive semi-definite. The determinant of this matrix is —2C oq 
attractive two-nucleon force where C < the conditions are satisfied provided that the 
three-body interaction coefficient D is not too large. We note that the positivity condition 
is spoiled more easily in the Hamiltonian limit where at — > 0. 



C. Instantaneous free pion action 

Before discussing lattice actions for chiral effective field theory, we first consider the lattice 
action for free pions with mass m n and purely instantaneous propagation, 

SmrjlTi) = 0^(^ + 3) ^ 7T/(n, 7lt)7Tf (w, Tl t ) - Ot t n I&> n *) 7r /(^ + 1 > n t) ■ ( 4 - 12 ) 

n,nt,I n,nt,I,l 
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The pion field 717 is labelled with isospin index /. Pion fields at different time steps rit 
and n' t are not coupled due to the omission of time derivatives. This generates instanta- 
neous propagation at each time step when computing one-pion exchange diagrams. It also 
eliminates unwanted pion couplings contributing to nucleon self-energy diagrams found in 
earlier work 65|. Though we call it a pion field, it is more accurate to regard ttj as an 



auxiliary field which is used to reproduce the one-pion exchange potential on the lattice. If 
for example we wish to c onsider low-energy physical pions within the framework of chiral 
effective field theory [161| , these scattering processes can be introduced perturbatively using 
external pion fields and additional auxiliary fields to reproduce the corresponding Feynman 
diagrams at each order. 



Following the notation in Ref. 



1421 ] . it is useful to define a rescaled pion field, n'j, 



7Tj(n,n t ) = v&ifa"t)> 



(4.13) 



where 



In terms of 77, 



q-x = a t (ml + 6). 



SW(tt/) = o ^ 7rj(n,nt)7Tj(n,n t ) - — 



n,n t ,I 

and in momentum space we have 



i.k 



1 aiv^ 1 

> COS ft; 

2 q, V 



The instantaneous pion correlation function at spatial separation n is 



(4.14) 



■K , I {n,n t )<K , I {n + i,n t ), (4.15) 

n,nt,I,l 



(4.16) 



(4.17) 



where 



Ar(fc) 



1 - — Ei cos h 



(4.18) 



D. Chiral effective field theory on the lattice 



We define some lattice derivative notation which will be useful later. There are various 
ways to introduce spatial derivatives of the pion field on the lattice. The simplest definition 
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for the gradient of tt'j is to define a forward-backward lattice derivative. For example we 
can write 



1 



d 1 n>(n) = -[n>(n + l)-ir' I (n-l)]. 



(4.19) 



This is the method used in Ref. 65J. One disadvantage is that it is a coarse derivative 
involving a separation distance of two lattice units. We can avoid this if we think of the 
pion lattice points as being shifted by —1/2 lattice unit from the nucleon lattice points in 
each of the three spatial directions. For each nucleon lattice point ra nuc i e on we associate a 
pion lattice point n pion , 



n 



pion 



1- 1- 1- 
^■nucleon 9^ 



(4.20) 



Then we have eight pion lattice points forming a cube centered at n nU cieon, 



n 



pion; "-pion 



^-pion + 1; ^-pion + 2, f-pion + 3 



^pion + 1 + 2, "ft-Dion + 2 + 3, 7^ D ion + 3 + 1, T^nion + 1 + 2 + 3. 



'pion 



pion 



''pion 



(4.21) 



For derivatives of the pion field we use the eight vertices of this unit cube on the lattice to 
define spatial derivatives. For each spatial direction I = 1, 2, 3 and any lattice function /(n) 
we define 



(4.22) 



i (-iy i+1 f{n+u), v = v 1 i + v 2 2 + v£. 

^1,^2,^3=0,1 

For double spatial derivatives of nucleon fields along direction I we use the simpler definition, 



^f(n) = f(n + l) + f(n-l)-2f(n). 



(4.23) 



At leading order in chiral effecti ve fi eld theory, the first partition function and transfer 



matrix operator considered in Ref. 



142] was 



Z 



LOi 



Tr 



(4.24) 



where 



M LOl =: exp \ -H iTCC a t - ^Ca^U'^n) * - ^Cj^J] \pf' a 



n 



+ 



8/X 



GsxS^fii - n 2 ) pistil) p a s£{n 2 ) \ : , 



(4.25) 



Si, 82, 1 ni,n,2 
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and 



GstSiin) = ^A Sl 7Tj(n, n t )A 52 7r^(0,nt)^) (no sum on I) 

= Yq E E (-ir Sl (-l)^(^(^+^-^^K(0,^)). (4.26) 

^1,^2,^3=0,1 ^,^2, 14=0,1 

This leading-order transfer matrix, labelled Mj J o 1 , has zero-range contact interactions anal- 
ogous to the pionless transfer matrix in Eq. (14. 6p . The 0( a 4 ) -improved action was used for 

A second leading-order partition function and transfer matrix was also considered, 



2lo 2 = Tr 



(Mi 



L0 2 



(4.27) 



where 



Mi 



L0 2 



: exp < -H{ ree a t 



at 
21} 



2 2 1 

+ ItSt E E G SiSa(«l -^)P5x,j("i)p&,j("2) > : • 
The momentum-dependent coefficient function f(q) has the form 



(4.28) 



f(q) = fo 1 ex P 



-b E (1 - cosgj) 



2=1,2,3 



(4.29) 



and the normalization factor f Q is determined by the condition 



/o-7^E 



exp 



/=1,2,3 



cos ; 



(4.30) 



The coefficient 6 was determined by fitting to reproduce the correct average effective range 
for the two S"-wave channels. For small q the function /(g) reduces to a Gaussian function, 



f(q) « / lex P (~^ 2 ) 



This Gaussian smearing of the contact interactions in Mlo 2 was 
clustering instabilities at lattice spacing a = (100 MeV) -1 142 



(4.31) 

bund to remove four-nucleon 
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E. Chiral effective field theory with auxiliary fields 



Let us define the auxiliary-field action 



n,nt,I 



In terms of auxiliary and pion fields, the partition function for LOi is 
Zl 0i = J Dtt'jDsDs! exp [S^ - S^} 

x Tr{M LOl (7r' I ,s,s/,L t -l) x ••• x M LOl «, s, s T , 0)} 



where 



M LOl (7r},s, s/,rii) =: exp 



-i/ fr eeat + y/-Ca t ^ s(n, n t )p at ' a (n) 



n,l 



n,S,I 



and Dtt'jDsDsi is the functional measure, 

<is(n, Tit) 



Dtt'jDsDsj = \[ 



n,n t 



271 



n 

J n,nt,I 



d7Tj(n, n t )dsi(n, n t ) 
2^ 



The instantaneous free pion action was already defined in Eq. (14.121) . 
For the LO2 action we have 

z ^ = x / Dn ' iDsDsi exp ~ 5 « oa ] 

x Tr{M L o 2 (7rj,s,Sj,L t - 1) x ■ ■ • x M L02 K, s, sj, 0)} 
The functional form of the transfer matrices are the same, 

M^i^i, s, S^Tlt) = M LOl (7r}, s, si, nt), 
but for LO2 the auxiliary-field action has the non-local form 

n,n',nt 

+ s i{n,n t )r 1 (n-n)s I (n \n t ), 



n,n' ,nt 
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where the inverse function / 1 is defined as 



-iq-(n—n') 



l 3 y f® 



(4.39) 



F. Next-to-leading-order interactions on the lattice 



The lattice studies in Ref. 144l . Il45| considered low-energy nucleon-nucleon scattering 
at momenta less than or equal to the pion mass, m n . On the lattice the ultraviolet cutoff 
momentum, A, equals it divided by the lattice spacing, a. As noted earlier, serious nu- 
merical difficulties appear at large A in Monte Carlo simulations of few- and many-nucleon 
systems. In attractive channels unphysical deeply-bound states appear at large A. In other 
channels short-range repulsion becomes prominent, producing destructive sign or complex 
phase oscillations. The severity of the problem scales exponentially with system size and 
strength of the repulsive interaction. 

In order to avoid these difficulties the approach advocated in Ref. 



144 



1451 ] was to set the 



cutoff momentum A as low as possible for describing physical momenta up to m n . In most of 
the published work so far the value chosen was A = 314 MeV ~ 2.3m n , corresponding with 
a = (100 MeV) -1 . This coarse lattice approach is similar in motivation to the continuum 
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3l| 



low-momentum renormalization group approach using V\ ow k C 

For nearly all \q\ < A the two-pion exchange potential can be expanded in powers of 
q 2 /{4ml), 



Ami 



Ami + 4' 



-Mq) 



3 Ami 

3 Ami 



(4.41) 



A (V N T L p EP ) 



384vr 2 /, 



2 fA 

7T _ 



Ami (8^1 - igl - 1) + ~q 2 {Ug\ - \lg\ - 2) + O ( ( ^ 



64tt 2 /. 4 



[(? ■ o'i) (q ■ ^2) - q 2 (0*1 ■ ^2 



1 + fe 



(4.42) 



This expansion fails to converge only for q near the cutoff scale A w 2.3m 7r , and so there is no 
practical advantage in 
Instead we simply use 



practical advantage in keeping the full non-local structure of V^lq P at this lattice spacing. 



v LO = vw + v OPEP , 



(4.43) 



30 



Vnlo = V LO + Ay' ' + V ( 2 ) , (4.44) 

where the terms in Eq. (I4.42p with up to two powers of q are absorbed in the definition of 
the coefficients for AV^ and V^ 2 \ 

Before describing the NLO lattice interactions in AV^ and V^ 2 \ we first define lattice 
current densities for total nucleon number, spin, isospin, and spin-isospin. Similar to the 
definition of A/ in Eq. ( 14.22}) . we use the eight vertices of a unit cube, 

v = v x \ + z/ 2 2 + ^ 3 3, (4.45) 

for ^1,^2,^3 = 0,1. Let v{—l) for I = 1,2,3 be the reflection of the / th -component of z7 
about the center of the cube, 

V{-i) = i!+(l-2vi)l. (4.46) 
The / th -component of the SU(4)-invariant current density is defined as 

nf > a (n) = \ E (- l ) Vl+1 <^ + »H))Hi(n + V). (4.47) 

^1,^2,^3=0,1 8,^=0,1 

Similarly for spin current density, 

</(^) = i E E (-ir + %(n + ^-l))[crs] l ,a,An + ^, (4.48) 

vi>V2,va=0>l i,j,i'=0,l 

isospin current density, 

</(^) = \ E E ^- l T l+l <^ + *H)) [n] jr aiAn + (4-49) 

^1,^2,^3=0,1 i,j,j'=0,l 

and spin-isospin current density, 

= £ E E (-l)^ 1 4,(^ + ^(-0)[^]^[r4,a, J ,(n + z7). (4.50) 

^1,^2,^3=0,1 ,j'=0,l 

In Ref. |144j the next-to-leading-order transfer matrices M NLOl and M NL o 2 were defined 
by adding the following nine local interactions to the leading-order transfer matrices Mj J o 1 
and Mlo 2 - The two corrections to the leading-order contact interactions are 

AV = \AC : ^p at ' a (n)p at - a (n) :, (4.51) 

ft 

AV P = l -AC P : J2 pf' a (n)pf' a (n) : • (4.52) 
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At next-to-leading order there are seven independent contact interactions with two deriva- 
tives. These are 

V q 2 = - l -C q 2 : ^Tp^^vjy^n) :, (4.53) 



2 

n,l 



V IW = —C^ : J^pf ' a (^fpf' a (n) :, (4.54) 
n,i,i 

V S 2 iq 2 = -^C 5 2 jg 2 : ^2ps' a (n)tfps' a (n) :, (4.55) 

n,S,l 

V S 2j2 >q 2 = ~C S 2j2 tq 2 : Y pi?(n)tfPsf(n) : > ( 4 - 56 ) 

n,S,I,l 

V iq . S )z = ^C [q . S )2 : Y Y A sPs' a ( H ) Yl A s'Ps'' a ( H ) : > ( 4 - 57 ) 

n S S' 

Vi*,(q-S)2 = Ifp^.sy : Y Y A sPsf(n) Y A S'pij(n) :, (4.58) 

n,I S S' 

V iiqxS) . h = - % -C {iqxS) . k : Y £ W K ' a (n)A l/p a s> a (n) + nf/(n)A,/> a (n)] : . (4.59) 

H,l,S,l' 

G. Model independence at fixed lattice spacing 

In effective field theory calculations model independence is often tested by checking sen- 
sitivity on the cutoff scale A. At a given order the difference between calculations for two 
different cutoff scales Ai and A 2 should be no larger than the omitted corrections at the next 
order. On the lattice this test is problematic since the lattice spacing cannot be changed by 
a large amou nt du e to computational constraints. Instead a different approach was intro- 



duced in Ref. 144J to test model independence at fixed lattice spacing which we summarize 
in the following. 

The notation ^/Q"/ A " is used to denote two-nucleon operators with the following prop- 
erties. yQ n /& n j s a sum Q f i oca i two-nucleon interactions that is an analytic function 
of momenta below the cutoff scale A and scales as n or more powers of momenta in the 
asymptotic low-momentum limit. The term "quasi-local" is used to describe V"^ ' A " since 
the interactions are short-ranged. At fixed lattice spacing we may consider two different 
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lowest-order actions with interactions of the form 

V LOl = V} 0) + V OPEP + vf /A \ (4.60) 

V L02 = V} 0) + V OPEP + V 2 Q2/A2 , (4.61) 

where V^ 1 ^ and vf lh? are different quasi-local operators with at least two powers of 
momenta. Since the leading-order interactions are iterated non-perturbatively the contact 
terms and in general have different coefficients. However low-energy physical 
observables should agree up to differences the same size as the omitted contributions at 
next-to-leading-order. 

Similarly at next-to-leading order we may consider two different actions of the form 

V NLOl = V LOl + AV; {0) + Vf } + V ± Qi/A \ (4.62) 
Vkx) 2 = V LQ2 + A\/ 2 (0) + Vf } + vf /A \ (4.63) 

where Vf 4/A4 and V^ 1 ^ are different quasi-local operators with at least four powers of 
momenta. Low-energy physical observables should again agree up to differences the same 
size as the omitted contributions at the next order. 

This technique provides a method for testing model independence of the low-energy lattice 
effective theory without changing the lattice spacing. In principle however it is good to 
check model independence in multiple ways, including different variations for \/ Q "/ A ™ as well 
as changing the lattice spacing as much as allowed by computational constraints. 



V. TWO-PARTICLE SCATTERING ON THE LATTICE 



A. Cubic rotation group 

Lattice regularization reduces the SO (3) rotational symmetry of continuous space to the 
cubic rotational group SO (3, Z). This group is also known as the proper octahedral group 
and abbreviated as O. This lack of exact rotational symmetry complicates the extraction 
of partial wave amplitudes. SO(3, Z) consists of 24 group elements generated by products 
of 7r/2 rotations about the x, y, z axes. Since SO(3,Z) is discrete, angular momentum 
operators J x , J y , J z cannot be defined in the usual sense. Let R z (vr/2) be the group 
element for a n/2 rotation about the z axis. The SO (3) relation 

(5.1) 

33 



R s (tt/2) = exp 



2 



can be used to define J z . The eigenvalues of J z are integers specified modulo 4. J x and J y 
may be defined in the same way using (vr/2) and R y (vr/2). 

There are five irreducible representations of the cubic rotational group. These are usually 
called Ai, T 1; E, T 2 , and A 2 . Some of their properties and examples using low-order 
spherical harmonics Yl,l z (#, 0) are listed in Table [H The 2 J + 1 elements of the total 
angular momentum J representation of SO (3) break up into smaller pieces associ ated with 
the five irreducible representations. Examples for J < 5 are shown in Table U H- In 

TABLE I: Irreducible SO(3, Z) representations. 



Representation 


Jz 


Example 


A x 


0mod4 


^0,0 


Ti 


0,1, 3 mod 4 


{^,0,^,1,^1,-1} 


E 


0,2 mod 4 


IV Y 2 ,-2+Y 2 , 2 \ 


T 2 


1,2, 3 mod 4 


( V Y 2 ,-2-Y 2 ,2 v 1 
[Y 2 ,l, ^ ,^2,-1) 


A 2 


2 mod 4 


>3,2->3,-2 
V2 



TABLE II: SO(3, Z) decompositions for J < 5. 



SO(3) 


SO(3,Z) 


J = 


Ai 


J = 1 


T X 


J = 2 


E@T 2 


J = 3 


Ti © T 2 © A 2 


J = 4 


Ax@Tx®E@T 2 


J = 5 


Tx®Tx®E®T 2 



lattice QCD these irreducible representations have been used to classi fy glueba ll states [163 ] 
as well as predict the spectrum and properties of baryon resonances jl64 . 165 ] . 



B. Liischer's finite volume formula 



Liischer's finite volume formula 1661 . 11671 . Il68| ] relates the energy levels of two-body states 
in a finite volume cubic box with periodic boundaries to the infinite volume scattering matrix. 
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FIG. 2: Sum of bubble diagrams contributing to two-particle scattering. 



Recently Liischer's method has been studied and exten ded in a number of different ways. 



Several investigations have looked at asymmetric boxes 



169, 



170], while another considered 



small volumes where the lattice length L is sma 
have also been studies of moving frames 
windings around the periodic boundary 
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ler than the scattering len 



172] , Yukawa interactions jl 731 ] , pion-exchange 



gth 



15] . There 



1741 ]. modifications at nonzero lattice spacing 175], 
and tech niqu es to distinguish shallow bound states from scattering states using Levinson's 
theorem 176] . Several recent studies d erived finite volume formulas for systems of n bosons 
with short-range interactions 1771 . Il78 ] . 

Liischer's method can be summarized as follows. We consider one up-spin and one down- 
spin in a periodic cube of length L. The two-particle energy levels in the center-of-mass 



frame are related to the S"-wave phase shift, 



p cot 6 (p) = —S (77) , 



V = \7T 



Lp 
2^ 



(5.2) 



where S(rj) is the three-dimensional zeta function, 

6(A 2 - n 2 ) 



s(v) 



lim 

A— >oo 



AttA 



rr 



V 



(5.3) 



The S-wave effective range expansion gives another expression for the left-hand side of 
Eq. (JE2D, 



p cot 5q (p) ~ — - 



Oscatt 



2 r oP + 



In terms of 77, the energy of the two-particle scattering state is 



P 2 _ 7] /27T 
-frpole — — — I ~ 



m 



m 



(5.4) 



(5.5) 



For the case of zero-range interactions, the location of the two-particle scattering pole is 
calculated by summing the bubble diagrams shown in Fig. [2j The relation between C and 



35 



£ P oie is [128| 



li^i Y, , „ — 7^77: (5-6) 



Ca t l^co L 3 _ ^ e -E polc a t _ i + 2a t u{2irk/L) - a 2 u 2 (2nk/L) ' 



k integer 

where 



if) = — Y] (1 - cos p;) (5.7) 

rn 



m 

1=1,2,3 



for the standard lattice action. In this manner the coefficient C can be tuned to produce 
the desired scattering length a scatt at infinite volume. Higher-order scattering parameters 
can also be extracted in this way. However for zero-range interactions the characteristic 
scale of these higher-order parameters is the lattice spacing, and so higher-order scattering 
corrections are the same size as lattice discretization errors produced by broken Galilean 
invariance and other lattice effects. 



C. Spherical wall method 

While Liischer's method is very useful at low momenta, it is not so useful for determining 
phase shifts on the lattice at higher energies and higher orbital angular momenta. Fur- 
thermore spin-orbit coupling and partial-wave mixing are difficult to measure accurately 
using Liischer's method due to multiple-scattering artifact s pro duced by the periodic cubic 



boundary. A more robust approach was proposed in Ref. 179j to measure phase shifts for 



nonrelativistic point particles on the lattice using a spherical wall boundary. Similar tech 



niques have long been used in nuclear physics (see for example Problem 5-7 in Ref. [18Q |) 
dating back to early work on i?-matrix methods 18 1| . We summarize the method as follows. 

A hard spherical wall boundary is imposed on the relative separation between the two 
particles at some chosen radius -R wa ii. This boundary condition removes copies of the 
interactions produced by the periodic lattice. Viewed in the center-of-mass frame we solve 
the Schrodinger equation for spherical standing waves which vanish at r = -R wa n as indicated 
in Fig. [3l 

When the combined intrinsic spin of the two interacting particles is zero there is no mixing 
between partial waves. At values of r beyond the range of the interaction, the spherical 
standing wave can be decomposed as a superposition of products of spherical harmonics and 
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FIG. 3: Spherical wall imposed in the center-of-mass frame. 



spherical Bessel functions. Explicitly we have 



[cos8 L ■ j L (kr) - sm5 L ■ y L {kr)] Y LiLz 



where the center-of-mass energy of the spherical wave is 

k 2 k 2 
E = 2 — = — , 



(5.8) 



2m 



m 



(5.9) 



and the phase shift for partial wave L is 5l- We can determine k from the energy E of the 
standing wave, and the phase shift 5l is calculated by setting the wavefunction in Eq. (15. 8p 
equal to zero at the wall boundary, 



cos5 L • j L {kR^\\) = sin6 L ■ y L {kR w& \\) 
j L (kR wsl i) 



S, 



tan 



(5.10) 
(5.11) 



_y L (kR wa ii) 

On the lattice there is some ambiguity on the value of R wa \\ since the components of r must 
be integer multiples of the lattice spacing. The ambiguity is resolved by fine-tuning the value 
of -R W aii for each standing wave so that 5l equals zero when the particles are non-interacting. 

When the combined intrinsic spin of the two interacting particles is nonzero, spin-orbit 
coupling generates mixing between partial waves. For nucleons the interesting case is S — 1 
where there is mixing between L = J — 1 and L = J + 1. We discuss this case here using 
the two- component notation, 

r Rj-i(r) 
Rj+i(r) 



(5.12) 
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for the radial part of the wavefunction. Since we are considering a two-channel system, 
there are two independent standing wave solutions of the form 

A T j_ x sin (k L r - ^-n + Aj_ x ) 
AS +1 sm(^r-^i7r + AS +1 ) 



V 1 oc 



kh 



(5.13) 



at energy E 1 = (k 1 ) 2 jvn and 



k n r 



(5.14) 



A 1 /^ sin (k H r - ^±tt + A 1 /^) 
_ A'/ +1 sin (k"r - ^ivr + A'/ +1 ) 

at .E^ = (k II ) 2 /m. These can be used to derive the phase shifts <5j_i and 5j + i and mixing 
angle ej using [179 ] 



tan (-Aj_ x + 5j-i) tan (— Aj +1 + Sj+i) = tan 2 £j, 

tan (-Ay_ x + (Jj.J tan (-Af/ +1 + 5 J+1 ) = tan 2 £j , 

1 1 s in (-Aj+i + 



cos 



Aj^tanej — — ^4j +1 - 

sin (-Af/ +1 + <f J+1 ) 



(-A^ + ^0 



A^i tanej 



-.4 



cos 



(5.15) 
(5.16) 
(5.17) 

(5.18) 



The phase shifts and mixing angle in Eq. (I5.15P and (I5.17P are at momentum k 1 while the 
phase shifts and mixing angle in Eq. (I5.16P and (15.181) are at momentum k 11 . Nearly equal 
pairs k 1 « k 11 are used in solving the coupled constraints Eq. (I5.15p - fl5.18p . In practice 
this amounts to considering the (n + l) st -radial excitation of L = J — 1 together with the 
n th -radial excitation of L — J + 1. Then we use 



tan (-A^_! + Sj^k 1 )) tan {-A T J+1 + <W^)) = tan 2 [ej^k 1 
tan (-Ay_ x + 5j„ 1 (k 1 )) tan (-A^ + <Jj+i(fc J )) ~ tan2 M fc/ 

1 J J+1 cos(-A',_ 1 +* J _ 1 (F)) 
for the phase shifts and mixing angle at k = k 1 , and 

tan (-Aj_ 1 + 5j-i(k n )) tan (— Aj +1 + <W*; J/ )) « tan 2 [ej(/c /7 )] 

tan (-Aj / _ 1 + ^_i(A; J/ )) tan (-Ay +1 + 5 J+1 (A; // )) = tan 2 [^(A; 77 )] 

sin (-Af/ +1 + 5 J+ i(^ 7 )) 



^ tan [£,(£;")] = -A*f +1 



cos(-Ay_ 1 + 5 J _ 1 (^))' 



(5.19) 
(5.20) 
(5.21) 

(5.22) 
(5.23) 
(5.24) 



for the phase shifts and mixing angle at k = k 



li 
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J z mod 4 



FIG. 4: Energy levels for S = 0, I = 1 using lattice actions LOi and LO2 and a spherical wall at 
radius i? wa ii = 10 + e lattice units 144] ] . The s olid line indicates the exact energy levels which 



reproduce data from the partial wave analysis of 



1821 ] . 



D. Scattering at NLO in chiral effective field theory 

Lattice phase shifts and mixing angles at leading order and next-to-leading order were 

100 MeV)-\ 



144| using the spherical wall method at lattice spacings 



a 



calculated in Ref. 

a t = (70 MeV) -1 . We summarize the results here. Fig. H] shows energy levels for spin 
S = and isospin 1 = 1 using lattice actions LOi and LO2. The spherical wall is at radius 
-Rwaii = 10 + e lattice units where e is a small positive number. The e notation makes explicit 
that |r| = 10 lattice units is inside the spherical wall but all lattice sites with \r\ > 10 lattice 
units lie outside. The solid lines indi cate the exact energy levels which would reproduce 



data from the partial wave analysis in 1821 ]. The energy levels for the standard action LOi 
are 10% to 15% too low for the 1 S'o states, while the improved action LO2 is correct to a 
few of percent for all 1 S'o states. Deviations for higher partial waves are smaller than one 
percent for both LOi and LO2. 

The energy levels for spin S = 0, isospin 1 = 0, and -R wa ii = 10 + e lattice units are shown 



39 



L0 1; S = 0,7=0 



LO 2 ,S = 0,/ = 



> 



:]::::::&::::: 



i]:::::: 



> 



o 

+ 
□ 



J z mod 4 



A 

E 
T 2 A 
A 2 X 
PWA93 — 



:]::::::&::::: 



W'.'.'M: 



J z mod 4 



FIG. 5: Energy levels for S = 0, I = using lattice actions LOi and LO2 and a spherical wall at 
radius i? wa ii = 10 + e lattice units 1441 ] . The s olid line indicates the exact energy levels which 



reproduce data from the partial wave analysis of 



1821 ] . 



in Fig. In this case LOi is better for the l P\ states and is within one percent of the exact 
values. The LO2 energy levels are further away, though still within a few percent for the 
1 Pi states. 

In Ref. 144j the nine unknown operator coefficients at next-to-leading order were deter- 



mined by matching three S'-wave scattering data points, four P-wave scattering data points, 
as well as the deuteron binding energy and quadrupole moment. Each of the next-to- 
leading-order corrections were computed perturbatively. The S'-wave phase shifts for LOi 
and NLO 1 versus center-of-mass momentum pcm are shown in Fig. El and the S-wave phase 



shifts for LO2 and NLO 2 are shown in Fig. 



71 The NLOi and NLO2 results are both in good 
agreement with partial wave results from 1821 ] . Systematic errors can be seen at momenta 
greater than about 80 MeV and are larger for NLOi. But in both cases the deviations are 
at larger momenta and consistent with higher-order e ffect s. 



P-wave phase shifts are shown in Fig. [S] and M 



i3- 



In this case the phase shifts 
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FIG. 6: 5-wave phase shifts versus center-of-mass momentum for LOi and NLOi 
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FIG. 7: S- wave phase shifts versus center-of-mass momentum for LO2 and NLO2 144]. 



are already close for LOi and quite accurate for NLOi. This suggests that only a small 
correction is needed on top of P-wave interactions produced by one-pion exchange. The 
results for LO2 and NLO2 are not quite as good. The Gaussian smearing introduced in LO2 
produces attractive forces in each P-wave channel that must be cancelled by next-to-leading- 
order corrections. However the residual deviations in the NLO2 results appear consistent 
with effects that can be cancelled by higher-order terms. 



The mixing parameter e\ for J = 1 is shown in Fig. [TO] 144]. The mixing angle is defined 
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FIG. 8: P-wave phase shifts versus center-of-mass momentum for LOi and NLOi 144 ] 



according to the Stapp parameterization 



1831 ] . Results for LOi and NLOi are on the left, 



and results for LO2 and NLO2 are on the right. The pairs of points connected by dotted 
lines indicate pairs of solutions at k = k 1 and k = k 11 for the coupled 3 S\- 3 Di channels. For 
LOi we note that E\ has the wrong sign. This suggests that the mixing angle may be more 
sensitive to lattice discretization errors than other scattering parameters. However for both 
NLOi and NLO2 results the remaining deviations appear consistent with effects produced 
by higher-order interactions. 
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FIG. 9: P-wave phase shifts versus center-of-mass momentum for LO2 and NLO2 144 ] 



VI. MONTE CARLO ALGORITHMS 



A. Worldline methods 



In bosonic systems or few-body systems where the problem of fermion sign cancellation is 
not severe, lattice simulations can be performed by directly sampling particle worldline con- 
figurations. We sketch an example of a lattice worldline configuration for two-component 
fermions in one spatial dimension in Fig. (TTTT) . This technique was used in the simulation 
of the triton using pionless effective field theory 149| . A number of efficient cluster algo- 
rithms have been developed for condensed matt er applications to gene r ate new worldline 



configurations based on loop and worm updates [184 



185 



186, 



187 



188 



1891 ]. 
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FIG. 10: £\ mixing angle for LOi and NLOi on the left, LO2 and NLO2 on the right [1441 ] 
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FIG. 11: Example of a worldline configuration for two-component fermions. 



While there are techniques which address the sign problem in certain cases 186], there 
is no general method known for eliminating sign oscillations in fermionic systems due to 
identical particle permutations. For Monte Carlo simulations extending over Euclidean 
time t, the sign of the configuration, sgn(C), averaged over all configurations C scales as 



(sgn(C)) ~ exp [(£ b ° SOniC - £fermiorn« 



(6.1) 



where ^ ermiomc i s the physical ground state energy and EQ OBomc is the fictitious ground 
state energy for bosons with the same interactions. The severity of the problem scales 
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exponentially with the size of the system and inverse temperature. In nuclear physics the 
same issue arises in continuous-space worldline methods such as Green's Function Monte 
Carlo and auxiliary- field diffusion Monte Carlo. In each case some supplementary condition 



is used to fix fermion nodal boundaries or constrain the domain of path integration 



190, 



191). 



B. Determinantal diagrammatic methods 



Determinantal diagrammatic Monte Carlo was used in 124j, Il25| to study two-component 
fermions in the unitarity limit near the critical point. This method is structurally similar 
to loop and worm updates of worldlines, however each configuration involves a complete 
summation of diagrams for a given set of vertices in Euclidean space. We discuss the 
method briefly here. 

Let be the free-particle propagator in Euclidean space. We note that G^ is real- 
valued. We define a set of n vertex locations 

Sn = m,t j )} j=h ... jn , (6.2) 

where fj is the spatial location and tj is the Euclidean time for the j th vertex. We also 
define a matrix of vertex-to- vertex propagators A[<S n ], where 

A i3 [S n } = G (0) (r1 - fj, U - tj). (6.3) 

As an example we choose a set of five points <Ss, and in Fig. (fl2|) we draw a Feynman 
diagram with vertices located at the coordinates of S$ . The propagators for the down spins 
in Fig. (FI21) give one term in the expansion of det A[«S 5 ], 

det A[S 5 ] = ■■■ + A U A 45 A 53 A 32 A 21 + ■■■ . (6.4) 

The same is true for the up spins in Fig. (1121) . 

det A[S 5 ] = A 25 A 53 A 32 A 14 A 41 + ■■■ , (6.5) 

and the determinant expansion shows that there is a relative minus sign between the up 
and down contributions. From this example it is clear that the total contribution of all 
Feynman diagrams with vertices given by S n is 

dP[S n ] = (-Ca t ) n {det A[S n ]} 2 . (6.6) 
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FIG. 12: One diagram contributing to the sum of diagrams with vertices located at the coordinates 
of S 5 . 



We note that dP[S n ] is positive definite when the interaction is attractive, C < 0. Conver- 
gence of the series in powers of C is guaranteed by finiteness of the Grassmann path integral 
at finite volume. 

In order to compute the full path integral 



2 = E / dp ^ 



(6.7) 



the sampling of vertex configurations can be generated using a worm algorithm that produces 
closed loop diagrams such as Fig. (fT2l) as well as single worm diagrams such as the example 
shown in Fig. ([TBI . In this diagram pairs of fermion lines are created at vertex 3 and 
annihilated at vertex 2. The sum of all diagrams of the type shown in Fig. f|T3|) can be 
written in terms of the derivative of det A[iS 5 ] with respect to A 32 , 

2 







OA 



32 



det A[S 5 



(6.8) 



From this we see that the contribution from worm diagrams is also positive. These diagrams 
are used to calculate the expectation value of the pair correlation function, 



(c l (r2,t2)c^(r 2 ,t 2 )c*(f 3 ,t 3 )cl(f 3 ,t 3 )) . 



(6.9) 



Further details of the worm updatin g algorithm and determinantal diagrammatic Monte 



Carlo can be found in 



124 



125, 



192]. 
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4 



FIG. 13: Single worm diagram with pairs of fermion lines created at vertex 3 and annihilated at 
vertex 2. 



C. Projection Monte Carlo with auxiliary field 



Projection Monte Carlo was usee 



fermions at unitarity [137 



139 



141 



;o co mpute the ground state energy of two-components 



143| . It was also used to study light nuclei and dilute 



neutrons in chiral effective field theory 1421 . 1 1451 ] . We briefly describe the method here 
using first the example of zero-range attractive two-component fermions. 

Let E^ N be the ground state for the interacting system of N up spins and iV down 



spins. Let 



0,free 
N,N 



be the normalized Slater-determinant ground state for a non-interacting 



system of N up spins and N down spins. We use the auxiliary-field transfer matrix defined 
in Eq. (I3.44p to construct the Euclidean-time projection amplitude 



J N,N 



w=nf/ 

n,n t 



d A s(n, n t ) 



0,free 
N,N 



M A (s,L t -l] 



M A (s, 0) 



0,free 
N,N 



(6.10) 



where t = L t a t . We define i?jv,iv(£) as the transient energy measured at time t, 

at) 



&N,N{t) - — In — -— 

oc t 6 NN (t) 



(6.11) 



So long as the overlap between 



0,free 
N,N 



and the ground state of the interacting system is 



nonzero, the ground state E% N is given by the limit 



N,N 



lim E N>N (t). 

t— >oo 



(6.12) 



As a result of normal ordering, MA{s,n t ) consists of single-particle operators interacting 
with the background auxiliary field and contains no direct interactions between particles. 
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Therefore we can write 



0,free 
N,N 



M A (s,L t 



M A (s,0) 



0,free 
N,N 



[det M A (s,t)f 



(6.13) 



where 

[M A (M)] W = (p k >\M A (s,L t - 1) M A (s,0) \p k ) , (6.14) 

for matrix indices k, k' — 1, • • • , N. \pk) , [py) are single-particle momentum states compris- 
ing the Slater-determinant state ^^T^. The single-particle interactions in M A (s,n t ) are 
the same for both up and down spins. Since the matrix is real- valued, the square of the de- 
terminant is nonnegative and there is no problem with sign oscillations. New configurations 
are accepted or rejected according to the probability weight 



dP(s) = [det M A (s, L t a t )} 2 d A s. 



(6.15) 



We note that \Ni A (s ) t)] k , k is only an N x N matrix. This is considerably smaller than 
matrices encountered in most other determinantal methods and contributes to the relative 
efficiency of projection Monte Carlo. For the case when the auxiliary field s is continuous, 
new configur ations can be generated using a non-local updating algorithm called hybrid 
Monte Carlo [1931 . Il94i . [195J . This scheme is widely used in lattice QCD simulations. 



D. Hybrid Monte Carlo 

We describe the hybrid Monte Carlo algorithm in general terms for probability weight 

P{s) oc exp [-V(s)} , (6.16) 

which depends on the lattice field s(n, n t ) and some function V(s) which may be a non-local 
function of s. The method proposes new configurations by means of molecular dynamics 
trajectories for 

n,n t 

where p(n,nt) is the conjugate momentum for s(n,n t ). The steps of the algorithm are as 
follows. 

Step 1: Select an arbitrary initial configuration s°. 
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Step 2: Select a configuration p° according to the Gaussian random distribution 



P [p°(n, n t )] oc exp <j -- [p°(n, n t 



Step 3: For each n, n t let 



~o / -> \ o / -» \ st °P 
p (n,n t ) = p — 



for some small positive £ s t op - 
Step 4: For steps z = 0, 1, iV s tep — 1, let 



s t+1 {n, n t ) = s\n,n t ) + e 8te pP\fi, n t ) 



p l+1 (n, n t ) = p l (n, n t ) - e step 



ds(n, nt) 



for each n, nt- 



Step 5: For each ft, n t let 



p Nstcp (n, m) = p NBtcp (n, n t ) + 



Step 6: Select a random number r G [0, 1). If 



dV(s) 
_ds(n, n t )J s=s iv stcp 



< exp [-H(s N «»,p rf '**>) + #(s°,p )] 



(6.18) 



(6.19) 



(6.20) 
(6.21) 



(6.22) 



(6.23) 



then set s° = s Nstcp . Otherwise leave s° as is. In either case go back to Step 2 to 
start a new trajectory. 

E. Grand canonical simulations with auxiliary field 



In Eq. ( 13.451) we introduced the partition function for zero-range attractive two- 
component fermions at chemical potential p, 



2(/*)=n j d 

n,n t 



l A s{n,n t 



Tr {M A (s, Lt-l,fi) M A (s, 0, //)} , (6.24) 



with auxiliary-field transfer matrix 



M A (s, n t , n) = M A (s, nt) exp <^ fxa t ^ p atfl ( 

I n 



n 



(6.25) 
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Let \n) be the quantum state with one particle at lattice site n and no other particles. As 
in Eq. (I6.14p . we define the one-particle matrix amplitudes 



[Ma(s, t)\ wn = (n'| M A (s, L t -1) M A (s, 0) |n) . 

However in this case the matrix [My^s, t)] H , H has dimensions L 3 x L 3 . 
The trace over states in Eq. (16.241) can now be written as 

Tr {M A (s, L t — 1, n) M A (s, 0, /i)} = {det [1 + e» at M A (s, L t a t )]f 



(6.26) 



(6.27) 



New configurations for s can be updated locally using the Metropolis algorithm. This 



method has been used in lattice calcu 
fermions near unitarity 



132 



133 



138 



ations to study the thermodynamics of two-component 



146 



147j | and, more generally, the attractive 



iubbard 



model and repulsive Hubbard model near half- filling in various spat ial dimensions 
A review of numerical aspects of this method can be found in 196 ]. 



59 



136]. 



F. Pseudofermion methods 



The same grand canonical partition function Z(/j,) in Eq. (16.241) can be evaluated in the 
Grassmann path integral formulation with auxiliary fields, 



z {v) = II J d A s(n,n t ) J DcDc* exp[-S A (c,c*,s, /i)] 



(6.28) 



where 



S A {c,c*,s,fi) = S A {e^c,c\s) + Yl K 1 ~ e^) c*(n,n t )c t (n,n t + 1)] , (6.29) 

n,n t ,i=hl 



and 



S A (c, c*,s) = S bee (c, c*) - ^ A [ s (™> n t)] P(™> n t)- 



(6.30) 



n,n t 



We note that S A (c,c* , s, /i) is a bilinear form coupling c and c* with a block-diagonal 
spin structure which is the same for up and down spins, 



S A (c,C*,S,fi) = J^XM ( # *' 71 *) [ Sa ( S > V)]n,n t ;H>,n> t c t(™'X) 
n,nt n' ,n' t 



(6.31) 
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Therefore the integration over Grassmann variables gives the square of the determinant of 



n,nt 



[n, n t 



{det S j4 (s,/i)} 2 . 



(6.32) 



This result can also be written as a path integral over a complex bosonic field (f)(n, n t ), 

d(j) ieal (n, n t )d(fii mag (n, n t 



z(n) = If 



n,n t 



d A s(n, n t )- 



2tt 



exp [-T A (<f), s, n)] 



(6.33) 



where 



T 



A 



n, n t 



S A (s,fi) S/ (s,fi) 



(6.34) 



The bosonic field <j)(n,n t ) is called a pseudofermion field. This technique was first im- 
plemented for fermions in lattice QCD 1971 ] . The non-local action in Eq. ( 16.341) can be 
updated using a non-local algorithm such as hybrid Monte Carlo. Typically an iterative 
sparse matrix solver is used such as the conjugate gradient method. 



Pseudofermion methods have 



fermions near unitarity 12 



j een used to study t he thermodynamics of two-component 



129 



130 



131 



132 



133| . For the case when an external field 



J is coupled to the difermion pair, 



22 [J*(n, n t )c\(n,n t )cl(n,n t ) + J{n,n t )ci(n,n t )c^(n,n t )\ , 



(6.35) 



n,n t 



the block structure of the Grassmann action is more complicated. However the analysis 



in Ref. 



60] shows that the path integral can still be written in terms of a positive-definite 



Pfaffian. Lattice simulations usi ng th is formalism were carried out using pseudofermion 
methods and hybrid Monte Carlo jl29| . 



G. Applications to low-energy nucleons 



The projection Monte Carlo method with 



energy nucleons in chiral effective field theory 1421 . 1 1441 . 1 1451 ]. A two-step approach was used 



aux iliary fie 



ds has been used to study low- 



where a pionless SU(4)-symmetric transfer matrix acts as an approximate and inexpensive 
low-energy filter at the beginning and end time steps. For time steps in the midsection, 
the full leading-order transfer matrix was used and next-to-leading-order operators were 
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operator insertion for 
expectation value 

<*| e w e | SU(4)/ fullLO (/ M1LO 5[/(4)/|^| e w e > 



2L, + L t . L t + L t . 

'o 'o '( 



L t +LJ2 

'o n 



FIG. 14: Schematic overview of the projection Monte Carlo calculation for nucleons in chiral 
effective field theory. 



evaluated perturbatively by insertion at the middle time step. A schematic overview of the 
transfer matrix calculation is shown in Fig. [TH 

The pionless SU(4)-symmetric transfer matrix is computationally inexpensive because 
the path integral in the SU(4) limit is strictly positive for a ny even number of nucleons with 
either spin-singlet or isospin-singlet quantum numbers 1601 ] . Although there is no positivity 
theorem for odd numbers of nucleons, sign oscillations are relatively mild in odd systems 
which are only one particle or one hole away from an even system with no sign oscillations. 



Some general results on positivity of t he path integra 



effective theory have been discussed in 



160 



198 



199 



and spectral inequalities in pionless 



200 ] 



201 



202], and the 



SU(4) symmetry arises naturally in the limit of large number of colors 
fact that both the spin-singlet and spin-triplet nucleon scattering lengths are u nusually large 
suggests that the physics of low-energy nucleons is close to the Wigner limit 2031 . 12041 ] . In 



Ref. 



2051 ] a general theorem on path integral positivity was derived for interactions governed 



by an SU(2iV)-invariant two-body potential whose Fourier transform is negative definite. It 
was also shown that as a consequence of the path integral positivity, the particle spectrum 
must satisfy a number of convexity lower bounds with respect to particle number. In 
Fig. (TI5T) we draw all SU(4) convexity bounds applied to the spectrum of light nuclei with 
up to 16 nucleons 205]. We note that each of the lower bound constraints are satisfied. 
While these results do not imply that Monte Carlo simulations of nucleons using chiral 
effective theory can be performed without sign or phase oscillations, they do suggest that 
the simulations are possible with only relatively mild cancellations given the approximate 
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nucleus 

n nn 3 H 5 He 6 He 7 Li ^elOBeUB 13 C 14 C 15 N 

'H 2 H 3 He 4 He 5 Li 6 Li 7 Be 8 Be 9 B 10 B l \C l2 C 13 N 14 N 15 p 16 
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(MeV) 



FIG. 15: Plot of the energy versus particle number for light nuclei with up to 16 nucleons. The 
line segments show the convexity lower bounds in the _SU(4) limit which hold for any two-body 



potential whose Fourier transform is negative definite 



205[]. 



SU(4) symmetry and attractive interactions at low-energies. 

VII. SOME RECENT RESULTS 

A. Ground state energy at unitarity 

At zero temperature there are no dimensionful parameters in the unitarity limit other 
than particle density. For up spins and down spins in a given volume we denote 
the energy of the unitarity- limit ground state as N . For the same volume we call the 
energy of the free non-interacting ground state E^f™^ and define the dimensionless ratio 

The parameter £ is defined as the thermodynamic limit for the spin-unpolarized system, 

£ = lim £ N)N . (7.2) 
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Several experiments have measured £ using density profiles of 6 Li and 40 K ex pand ing upon 
relea se from a har moni c trap. Some recent measured values for £ are 0.51(4) 2061 ] . 0.46^qs 



2071 ]. and 0.32+}° |208(. There is some disagreement among these r ecen t measure ments as 



209, 



210 



21l| 



well as with larger values for £ were reported in earlier experiments 

There are a number of analytic calculations for £ using techniques such as BCS sad- 
dle point and variational approximations, Pade approximations, mean field theory, density 



functiona 



pansions 



the ory, exact renorma l ization group, dimensional e-expansion s, and large- N ex 



212 



213 



214 



215 



216 



217 . 



218 



219 



220 



221 



222 



223 . l224j . The values for £ 



range from 0.2 to 0.6. Fixed-node Green's fu nctio n Monte Carlo simulations for a per iodic 



cube find £ NN to be 0.44(1) for 5 < N < 21 |225| and 0.42(1) for l arger N [22$ [227|. A 



228], and a mean- field 



restricted path integral Monte Carlo calculation finds similar results 
projection lattice calculation yields 0.449(9) |l43| . 

There have also been simulations of two-component fermions on the lattice in the unitarity 
limit at non zero temperature. When data are extrapolated to zero temperature the results 



of 



124 



38, 



3 



1461 ] produce a value for £ similar to the fixed-node results. The same is true for 



1251 ] . t hough wit h significant error bars. The extrapolated zero temperature lattice 



results from 



130 



13JJ established a bound, 0.07 < £ < 0.42. 



R ecent lattice calculations in the grand canonical ensemble yield a value for £ = 0.261(12) 



132 



1331 ] . These calculations used lattice volumes of 4 3 , 6 3 , 8 3 , 10 3 and also probed the 



jeha vior at finite scattering length. In Fig. (fT6i) we show £ as a function of rj — k F a 



i„-i 



scatt 



1331 ] . The circles show the lattice results of [1331 ] , and the dotted line shows a quadra tic fi t 



through the points. The squares are fixed-node Green's function Monte Carlo results 



and the so lid li ne corresponds with results calculated using the epsilon expansion [229 1. 



226j, 



In Ref. |139l ] £n,n was calculated on the lattice using Euclidean time projection in small 
volumes where it was estimated that £ = 0.25(3). More recent results using a technique 
called the symmetric heavy-light ansatz found similar values for £n n at the same lattice 



volumes and estimated £ = 0.31(1) in the continuum and thermodynamic limits 230| . A 
very recent lattice calculation using Euclidean time projection with a bounded continuous 
auxi liary field used lattice volumes 4 3 , 5 3 , 6 3 , 7 3 , 8 3 and extrapolated to the continuum limit 



137l | . The results found were 



£ 5 , 5 = 0.292(12), 
£ 7i7 = 0.329(5). 



(7.3) 
(7.4) 
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FIG. 16: Plot of £ as a function of r\ = k p a scatt . The circles show the lattice results of 



squares are fixed- node Green's fun ction Monte Carlo results 
with epsilon expansion results 2291 ] . 



133]. The 



2261 ] . and the solid line corresponds 



In Fig. [T7] we show resu lts for £ 5)5 and £7^ at finite L and the corresponding continuum 



limit extrapolations 



137]] . For comparison we also show Hamiltonian latt ice results using 



23oj . 



These lattice 



the symmetric heavy-light ansatz in the lowest filling approximation 
calculations show close agreement with each other and disagreement with fixed-nod e Gr een's 
function Monte Carlo results for the same number of particles in a periodic cube 225 ] . 



B. Critical temperature at unitarity 

At unitarity the critical temperature T c can be written as a fraction of the Fermi e nergy . 
Experimentally T c /Ep has been measured using trapped 6 Li and found to be 0.27(2) 206 ]. 



However the interpretation of this result is made difficult by modifications caused 



potential and the problem of relating empirical and actual temperature scales 



j y the trap 



124 



231 



232 ] 



A number of approx imate theoretica l calculations s uggest a value for the critical temperature 



above 



71 



231 



2331 ] as well as below 234], 12351 . 12361 ] the Bose- Einstein condensation tempera- 



ture Tbec — 0.218-Ef. An epsilon expansion calculation around d = 2 yields Tr /Ep w 0.15 



218 

3 



237| . Omitting 



22J. A continuum-space 



while the epsilon expansion around d = 4 yields T c /Ep pa 0.25 
terms at 0(N~ 2 ), the large N expansion yields T c /Ep ~ 0.14 
restricted path integral Monte Carlo calculation found T c /Ep ~ 0.25 |228 |. 

Lattice simulations measuring long-range order in the pair correlation function find values 
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? IG. 17: Results for £5 5 and £7 7 at finite L and the corresponding continuum limit extrapolations 



1371 ] . For comparison we also show Hamiltonian lattice results using the symmetric heavy- light 



ansatz in the lowest filling approximation 



2301 ] . 



124 



125), and T C /E F 



TJE F < 0.14 [131|, T C /E F < 0.15(1) [143], TJE F = 0.152(7) 
0.183(12) 133[ | . The spread in values can likely be explained by lattice discretization 
errors, which are visible in Fig. ffTBl) showing the dependence of T C /E F on t> 1//3 , where v is 



124 



the lattice filling fraction 
12 3 . The point labelled A. Sewer et al. corresponds with 



1251 ] . The simulations were done with lattice sizes 6 3 , 8 3 , 



59(, while the points labelled T. 



A. Maier et al. corr espond with unpublished work which appears to be unavailable in print. 



The results of 



1291 ] are also consistent with a point along this line. 



While coherence measurements of the pair correlation function in 



147] indicate an upper 



bound on the critical temperature, T c f E F < 0.15(1), the calculatio n of the a verage energy 



138 



1471 ] . This data is 



has a peculiar structure at T/E F = 0.23(2) at lattice volumes 6 3 , 8 3 

shown in Fig. (JT9l) . The physical significance of this effect is presently unknown. Meanwhile 
lattice calculations of the pair correlation function using projectio n M onte Carlo find low- 
energy string-like excitations winding around the periodic lattice 1411 ] . These excitations 
may play some role in spoiling pair coherence at relatively low temperatures. 
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FIG. 18: The critical temperature T c /Ep versus v 1 ^ 3 , where v is the lattice filling fraction 
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1251 ] . 
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IG. 19: Plot of the average energy per particle in units of \Ep versus temperature in units of Ep 



138 



1471 ] 



C. Dilute neutron matter at NLO in chiral effective field theory 



In Ref . 



145J the ground state energy of dilute neutrons was calculated on the lattice at 



next-to-leading order in chiral effective field theory. The simulations used 8 and 12 neutrons 
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FIG. 20: Results for Eq nlo/^q 6 versus Fermi momentum fcp 



145]. 



in lattice volumes 5 3 , 6 3 , 7 3 at lattice spacings a = (100 MeV)" 1 , a t = (70 MeV) -1 . In 
Fig. [20] we show results for the ratio of the interacting ground state energy to non- interacting 
ground state energy, Eq^lo/Eq^, as a function of F ermi momentum kp. For comparison 



we show other results 



v 



rum the literature: I I' 1981 1238], APR 1998 |239j, CMPR v6 and 
240J, SP 2005 |241| . and GC 2007 [39j. There is good agreement between the different 
resul ts ne ar kp = 120 MeV, but there is some disagreement on the slope. The analysis in 
Ref. 145j shows that much of the P-wave contributions from different spin channels cancel 



numerically 



D. Comparison with other methods and future outlook 



At nonzero temperature there are unfortunately very few ab initio calculations that can be 
used to compare with results obtained using lattice effective field theory. We have already 
ment ioned a restricted path integral Monte Carlo calculation for cold atoms at unitarity 



228( | . However the size of systematic errors due to path restriction is difficult to estimate, 



and the final result for the critical temperature is in strong disagreement with each of the 
lattice results presented above. 

More comparisons can be made for calculations of low energy spectra. At present the 
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most accurate ab initio calculations of light nuclei binding energies for up to twelve nucleons 
have been obtained using Green's Function Monte Carlo. The overall accuracy of these 
calculations are at the 1 — 2% level. Current lattice calculations are not as accurate as 
this, but it is hoped that lattice simulations for light nuclei at N 3 LO in the future can reach 
comparable accuracies. At cutoff momentum A = 500 MeV, No-Core Shell Model calcula- 
tions using the NNLO chiral potential (plus N 3 LO ter ms fo r the two-nucleon potential) give 



binding energies for Li and Li at 3% accuracy 



52 



242l | . The first lattice effective field 
theory calculations at NNLO are currently in progress. Preliminary results for l attice spac- 
ing a = (100 MeV) -1 give an alpha binding energy accurate at the 5% level 2431 ] . However 
much further work remains in developing the lattice formalism at higher order and studying 
larger nuclei. 

For the ground state of cold atoms at unitarity and dilute neutron matter, the quality of 
lattice effective field theory calculations are competitive with or exceed other computational 
approaches. Here the relative success of lattice simulations over other methods is probably 
due to the nature of the ground state and the use of determinantal Monte Carlo. For cold 
atoms at unitarity and dilute neutron matter the competition between attractive binding 
forces and Fermi antisymmetry is somewhat evenly matched, resulting in a complicated 
superfluid ground state somewhere in between BCS and BEC. This makes it difficult to 
use techniques where nodal constraints must be approximately guessed and easier to rely on 
determinantal Monte Carlo methods which automatically incorporate Fermi antisymmetry. 

In addition to probing more nucleons and higher orders in effective field theory, future 
work must also probe simulations at larger lattice volumes. This includes both smaller 
physical lattice spacings as well as larger physical volumes. The transition from small 
lattice systems to large production runs should be possible as more experience and data 
is collected on the efficiency of various lattice algorithms and more computing resources 
are devoted to important calculations. It is probably unlikely that lattice effective field 
theory simulations can match the spatial lattice lengths L = 30 ~ 40 used in some large- 
scale lattice QCD simulations. For cold atom simulations at unitarity the most significant 
computational barrier with increasing system size is the increase in condition number for the 
auxiliary field matrices M^(s,t). Similar computational slowdown occurs in unquenched 
lattice QCD simulations at very small sea quark masses or at large chemical potential (in 
addition to the problem caused by complex phase oscillations). 



59 



For a single bound nucleus it is not necessary to probe volumes much larger than the size 
of the bound state since the finite volume errors are exponentially small. In this case it 
would be more useful to probe smaller physical lattice spacings, as much as constraints such 
as sign or complex phase oscillations will allow. For unbound nuclear systems the finite 
volume dependence of energy levels should be more interesting. This data can be used to 
probe nucleon-nucleus scattering or nucleus-nucleus scattering using Liischer's finite volume 
scattering method. 

It is difficult to predict the development of the field in the future, but one general hope is 
that collaborative efforts develop with researchers not directly involved in large-scale lattice 
calculations. One working model has already been pioneered in the lattice QCD community, 
where large numbers of gauge configurations are stored and shared for general use. A similar 
model may be useful for lattice effective field theory calculations for systems with significant 
general interest. 

VIII. SUMMARY 

In this article we have reviewed the current literature on lattice simulations for few- and 
many-body systems. We discussed methods which combine the theoretical framework of 
effective field theory with computational lattice methods. The lattice spacing serves as the 
ultraviolet cutoff for the low-energy effective theory, and all interactions are included up to 
some chosen order in power counting. By increasing the order, the accuracy at low energies 
can be systematically improved. One feature of this approach is the ability to study several 
different phenomena using exactly the same lattice action. Another feature of the lattice 
effective field theory approach is the close theoretical link with standard analytic tools used 
in effective field theory calculations. The approach also benefits from the computational 
flexibility provided by a number of efficient lattice algorithms. We have discussed many of 
these in this article. 

The idea of lattice simulations using effective field theory is relatively new, and this 
review article represents a snapshot of the current progress in the field. We have attempted 
to cover the relevant principles from effective field theory as well as different formalisms 
and algorithms used in recent lattice calculations. We have focused much attention on 
techniques which can be applied to both cold atoms and low-energy nuclear physics as well 
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as common methods used in work by different collaborations. 
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